Periodic orbits in scattering from elastic voids 
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Abstract. The scattering determinant for the scattering of waves from several obstacles is consid- 
ered in the case of elastic solids with voids. The multi-scattering determinant displays contributions 
from periodic ray-splitting orbits. A discussion of the weights of such orbits is presented. 



1. INTRODUCTION 

Eventhough the word scattering appears to imply transport, studies of Helmholtz scat- 
terers have shown effects of trapped periodic orbits flj. As there is nothing particular 
about the Helmholtz equation, similar effects are expected for wave equations for other 
media such as dielectric or elastodynamic ones. 

We shall discuss the relation between periodic trapped rays and the scattering determi- 
nant corresponding to a medium with several polarizations, each with their own velocity. 
The example to be treated is the case of elastic wave propagation in a solid punctured 
by a finite number of voids. These systems have the feature that a ray hitting a bound- 
ary can either reflect or refract. Particularly, ray splitting occurs when the polarization 
changes. This leads to a ray dynamics which no longer is unique, since - in general - a 
single polarized ray evolves into a tree of rays. A similar behaviour is observed in mi- 
crowave resonators with dielectrica: characteristically rays can either be reflected or be 
transmitted at the boundaries of the dielectrica. 



The Helmholtz equation 



2. SCALAR CASE 



{A + k^)}j/ = (1) 



describes the wave propagation of a scalar field i// of wave number in a homogeneous 
and isotropic medium. Furthermore, if obstacles are embedded in this background, the 
scalar field typically has to satisfy Dirichlet (i// = 0) or Neumann {d\ir/dn = 0) boundary 
conditions on the surfaces of the obstacles. For an exterior problem, i.e. a scattering 
problem, the geometry of the scatterers and the corresponding boundary conditions are 
usually specified and the typical goal is to calculate the scattering matrix S. In short, this 
matrix contains the information on how an incoming wave transforms to a superposition 
of outgoing scattering solutions. From the knowledge of the scattering matrix several 



interesting quantities can be calculated: cross-sections, resonances, phase shifts, time- 
delays etc. 

One way of finding the scattering matrix is via the so-called null-field method 10, 
I2I S> Bl 1^ El d. 113] • A given field on the boundary of one scatterer gives rise to 
secondary fields on the full set of boundaries, including the one at infinity. These fields 
are calculated via boundary integral identities. If a basis is chosen for each boundary, 
the initial and the secondary fields can be expanded in these bases. The matrices that 
describe their relationships can be used to construct the scattering matrix Li this 
way the application of the standard null-field method determines the S matrix as 

S = 1-/CM"^D, (2) 

where the C and D matrices connect the incoming, respectively, outgoing waves to the 
interior scattering boundaries and where the matrix M relates the waves at one interior 
scattering boundary to the waves at a different one. M itself can be written in terms of a 
transfer matrix A as 

M = l + A. (3) 
The multi-scattering expansion arises when 

= 1-A + A2 + ... (4) 

is inserted in Q ||2,SB|5|]. This signals that M has the role of an inverse multi- scattering 
matrix. In some situations the exact matrix elements of these infinite dimensional matri- 
ces, A, C, and D, can be derived. Of course these matrices depend crucially on the scat- 
tering problem in question, but the overall structure is the same. In general, the null-field 
method can be applied to solve numerous multi- scattering problems in various media. 

An explicit example is the scattering of a scalar wave, propagating in the two- 
dimensional plane, off a finite number of non-overlapping hard discs that have fixed 
positions and Dirichlet boundary conditions. For the discussion here only the ma- 
trix M is of interest. The fields at two different discs j and /, expanded in basis 
states exp{im^j) and exp(«m'0y) (in polar coordinates with m,m' = 0, ±1, ±2, ■ ■ ■ two- 
dimensional angular momentum quantum numbers) are connected via the (transfer) ker- 
nel Jlliliil [3 El El 



[Ayr , = (1 _ 5„04#^^^'^ '(kRw) . (5) 

Here Rjji is the center-to-center distance of the discs of radii Uj, and aj/^ is the angle 
to the center of cavity / in the coordinate system of cavity j. All this follows from the 
application of the above- sketched null-field method to the scalar-wave problem. Less 
familiar is the ray limit of the pertinent multi-scattering determinant DetM(^) Jfl 0]. 
As the various building blocks of the scattering matrix are known in analytical form 
it is possible to explicitly calculate the short-wave-length limit. The main result is that 
asymptotically 



ir{kLi,-Vpn/2) 

DetM(.).F^exp|-E£-^^^^^-^^.-. 




(6) 

z=l 



where p is a prime orbit (a periodic orbit that cannot be split into shorter ones), r is 
the number of its repeats, Lp is its length, Vp is its Maslov index, Jp is its monodromy 
matrix, and Hp is its number of bounces on the scatterers. The right-hand side of © is 
the semiclassical spectral determinant . The parameter z, which at the end is put 

equal to unity, serves in the formal expansion of F in z up to a sufficiently high number 
Nmax >np. A periodic orbit is defined here as a closed ray in phase space obeying the 
law of reflection at the scatterers. Our goal is to generalize these types of problems to 
two-dimensional elastodynamics. 



3. ELASTODYNAMICS 



In isotropic elasticity the wave equation in the frequency domain CO is 

/iAu + (A + /i) V(V ■ u) + pco^u = 0, 



(7) 



where u(x) is the displacement vector field in the body. A, /i are the material-dependent 
Lame coefficients and p is the density [ l^ilTJ. This wave equation admits two different 
polarizations: longitudinal L and transverse T with velocities 



cl 



and 



Ct 



(8) 



The longitudinal and transverse waves correspond to pressure, respectively, shear defor- 
mations (or to primary and secondary arriving pulses in seismology). This leads to the 
law of refraction for incoming plane waves 



Ct 



sin 

sin 6t 



(9) 



where 6l, Ot denote the angle of incidence or reflection of the longitudinal and trans- 
verse wave, respectively, measured with respect to the normal to the surface. 
The stress tensor in elasticity has the form 



Oij = A dkUk5ij+ll [diUj + djUi) . 
The boundary conditions considered here are/ree. Hence 



t(u) = C7(u) -11: 



A(v-u)l + Ai{(Vu) + (Vu)'^} 



n = 



(10) 



(11) 



for the displacement field at the boundary where n denotes the normal to the boundary 
and T indicates a transposition. The operator t refers to the traction. 




FIGURE 1. General zone of two-dimensional cavity scattering. Here the unit vectors n denote the 
outside normals to the surfaces (9, (/ = 1,2,3) of the cavities or to the surface d„ of the total scattering 
zone. Furthermore, an incoming plane wave and a scattered cylindrical wave are shown. 



4. SCATTERING DETERMINANT 



As mentioned, in this treatment the medium corresponds to a cylindrical solid made, for 
simplicity, of an isotropic and homogeneous elastic material. The scattering geometry 
consists of parallel cylindrical voids, which are perpendicular to the endcaps of the 
overall cylinder, see fig.H] If the fields are stimulated by in-phase line-sources parallel to 
the voids, this symmetry is respected and the problem reduces to one of two-dimensional 
elasticity, see Fig. 1, referred to as plane strain. This scattering problem generalizes the 
simpler scalar scattering off discs Jn two dimensions, mentioned in sect.|2l 

Similar to the scalar case J^llOllllI], the scattering determinant may be factorized 
into an incoherent single-scatterer part (in terms of the determinants over the single- 
scattering matrices S^^-'^((d), j = 1, 2, ■ ■ ■) and the genuine multi-scatterer part (in terms 
of the determinant of the inverse multi-scattering matrix M((o)) [9] 

Uclk ^ 'i Det[M (<!>)] 

When the relative positions of the cavities are changed, only the latter factor, composed 
from the cluster determinant Det M((o), changes: 

M = l + A, (13) 



(7—La'=L 



Oa(y'n,_i,[kaKjf)e J ^ 



The indices G {r, 0} are the labels of the two-dimensional polar coordinates and 
the indices a, a' G {L,T} refer to the two polarization states. The diagonal matrix 

T^P'^'^ may be interpreted as a translation matrix acting on the polarized scattering 

states, where ki = (o/ci and fcj- = (o/cj are the pertinent wave numbers JvlH]- As in 

(jSj), Rjji is the center-to-center distance of the circular cavities of radii aj and Oj^ is the 

angle to the center of cavity / in the coordinate system of cavity j. The single-cavity 
scattering matrices (enumerated by the cavity index j = 1,2,...) are separable in the 
(two-dimensional) angular momentum Z = 0, ±1, ±2, ■ ■ ■ due to the rotational symmetry. 
They have the general form: 

[s!^)^]=-[tWr^-[f^'] (14) 

in terms of the 2x2 traction matrices [tj^^''] which incorporate the free boundary 
conditions (fTTT) . Here the "type" Z G {-f, — refers to outgoing, incoming or regular 

scattering states and involves H^^\h^^^^ or 7/ Bessel functions of argument zl,t = «^L,r> 
respectively. Thus we have, e.g., for the outgoing case: 

I {I -Zt/2)Hi {zl)-Zl 'azL 'M^/ (zr)-zr '^zt ) 



(15) 

Note that the single-cavity scattering matrix connects different polarizations. For a full 
discussion, s ee [i9nl8i.il9i.i2Qi1 . The connection to the interior problem of a single disc is 
described in 11211]. 

As first shown for the scalar problem |6], the poles of the cluster determinant 
DetM((o) cancel - by construction - the poles of the single- scattering determinants. 
Likewise, the poles of DetM((0*) are canceled by the zeros of the single-scattering de- 
terminants. Thus all scattering resonances defined by the poles of the scattering determi- 
nant detS((o) can be found from the zeros of the cluster determinant DetM(a)) flITlll. 
As an example consider the resonances in fig. El see 1I20I] . of a two-cavity system made 
of polyethylene [18], a material with ct = 1950 m/s and ct = 540 m/s. Furthermore, it is 
assumed that the cavity radii are equal, i.e. ai = a2 = a, and that the inter-cavity separa- 
tion R, measured from the centers, is 6 times larger than a. Note that the regular spaced 
horizontal set of resonances in fig.|2lis placed below an irregular set. This is opposite to 
the scalar Helmholtz case for the same geometry where the regular spaced resonances 
are above the irregular ones Ll_4. 22^ 23, 24]. The regular resonances particular to the 
fundamental A 1 -representation are well described by the following condition l[l2ll 

0=l+exp{ikLL)/VA (16) 



with the length L = 4a and instability A = 5 + 2^/6 JH]. L corresponds to the length 
of the shortest periodic orbit moving in a symmetry-reduced domain spanned by the 
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FIGURE 2. Elastodynamic scattering resonances from lEflll for two cavities of common radius a (Ai- 
representation, center-to-center-separation 7? ~ 6a) in the complex longitudinal Helmholtz number-plane: 
kia — COa I ci. 



surface of the cavity and the center-of-mass of the two cavities. A is obtained from the 
product of ray matrices as the leading eigenvalue of the monodromy matrix Jll, [Isl] of 
the corresponding (geometric acoustic) ray system. This next raises the question about 
the effect of the remaining set of orbits. 



5. ORBITS IN TIME-DELAY 

For real frequencies the total scattering phase is given by the sum over the cluster 
phase ©c and the single cavity phases 0^: 

1 1 DetA/ftt)*)^ 

0(«) = -lndet5(w) and 0,.(w) = -In . (17) 

^ ' 2i ^ ' ^ ' 2i DetM(£t)) 

see (fT^ . Likewise, the derivative with respect to frequency d&/dco, the Wigner-Smith 
time delay, can be decomposed into single- scatterer and cluster contributions. The nu- 
merics of the cluster time-delay d&c/dco show fluctuations which are related to trapped 
periodic orbits in the scattering geometry, see fig. |3] where the results for two identical 
cavities (the same system as in fig.EJi are presented. Due to the symmetry of the system 
the cluster delay decomposes further into a sum over four irreducible representations of 
the symmetry group C2v [25], of which the representations Ai and B2 are shown. Some 
of these orbits are diffractive (see [9, 20] for more details) including segments of sur- 
face propagation of Rayleigh type, which are also important in earthquakes lE^ . For 
these proceedings we focus the discussion to purely non-diffractive contributions, called 
geometrical ray- splitting orbits. 
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FIGURE 3. Power spectrum of the time-delays derived from the cluster phase shift of two cavities 
of common radius a and center- to-center separation R = 6a. The thick and thin line denote the Ai- and 
B2-representation, respectively. The symbols Pi and Si label, in turn, periodic orbits with / legs of pressure 
or shear polarization. The circular arcs indicate Rayleigh-surface waves. 



6. EXPANDING THE CLUSTER DETERMINANT 



As in the scalar case Jg], a central point in the orbit-construction is the definition of the 
cluster determinant in terms of traces: 



F(z) =Det(l + zA) =exp 



z"Tr(-A)" 



(18) 



where z is again a formal expansion parameter which is put equal to one at the end. 
Equation (HEJl holds if A is trace-class, and its Taylor expansion is called the cumulant 
expansion [6, 1 ll 12, S- Trace-class operators (or matrices) are those, in general, non- 
Hermitian operators (matrices) of a separable Hilbert space which have an absolutely 
convergent trace in every orthonormal basis iEtLEsIi . Especially the determinant Det ( 1 + 
zA) exists and is an entire function of z, if A is trace-class. Presently the trace-class 
property of A has only been proved in detail in the two-dimensional scalar case [6, 11] 
and sketched for the three-dimensional case in iTulh . Nevertheless, we shall proceed as if 
this is true also in our elastodynamical case. This is supported by the following numerical 
evidences: (i) the sum over the moduli of the eigenvalues of the matrix A from (fT3l) is 
absolutely converging in agreement with the expected trace-class property of this matrix; 
(ii) the determinant (fT8l) converges to a finite result as the dimension of A is increased 
beyond a minimal number A^dim > 2 x (§(cL/cr) I^lI^i), see (id] and also 10,0]; and (iii) 
the resonances of the cumulant expansion truncated at fourth order z^ agree very well 



with the exact ones plotted in fig. El see fid]. 



7. RAY LIMIT AND ORBITS 

The expansion (fTSt indicates that the cluster determinant can be obtained from the 
knowledge of an increasing number of traces. Moreover, each trace TrA" is given by 
the sum over all exact periodic itineraries of topological length n, see [6]. In the saddle- 
point approximation the periodic itineraries become the periodic orbits of topological 
length n, which means that they bounce n times between the cavities f^. These orbits 
fulfill the laws of reflection and refraction and have phases corresponding to their time 
periods of revolution Tp. 

A periodic itinerary of topological length n corresponds to a cyclic product of n terms 
involving one operator of the type 

[t!^^r'-[tp] (19) 

(which corresponds to the T-matrix-part (times —1/2) of a single-cavity scattering 
matrix ST4\i ) followed by one translation operator T^^t^"'"' (the free propagation). The 

pruning rule that two successive scatterings must take place at different cavities is 
automatically built in (see the term (1 — 5yy) in the kernel A). The ray limit of the single 
cavity T-matrix gives unitary reflection coefficients similar to those of the scattering 
from an infinite half-plane [16, 21]. This leads to an overall amplitude ttp defined as 
a product over all reflection coefficients along the orbit. This amplitude describes the 
leakage from the orbit due to ray splitting. 

The calculation of the geometric amplitudes of the orbits requires more work. See 1E9I] 
for a general discussion with respect to the interior scalar case and L6J for the exterior 



counter part. Asymptotic wave theory indicates 1130113 1LI32L 13311 that for open trajectories 
in two dimensions the amplitude scales as (fci?)~^/^ where k is the wave number in 
question and R is the radius of curvature of the wave front at the observer. This radius is 
studied in e.g. geometric optics. It it is possible to keep track of its evolution in the free- 
propagation period between the scatterers and at impacts including possible refractions 
with the help of suitable ray matrices [9] . Indeed, for our problem it can be shown that 
all those open segments that have fixed end points, but intermediate points (variables) 
determined by saddle-point integrations, have such an amplitude evolution. This comes 
about by calculating the accompanying sparse Hessian of this restricted integration. 

For a full saddle-point integration over all variables, in other words for a periodic 
orbit p, the amplitude turns out to be expressible as yet another sparse Hessian that can 
be expanded into Hessians of the type of the previously considered open pieces; see [6] 
for the scalar case. The use of the previous information then allows the full calculation 
with the amplitude evolving as 

^ - P ^n,^ (20) 



|Det(l-J;,)|V2 



where Jp is the product of the ray matrices and and Up is the product of the reflection 
coefficients, calculated along the orbit p with its np bounces. This form is precisely part 



of the conventional semiclassical density of states However, the formal 

parameter z is also present and can be seen as a counting and ordering parameter of the 
various orbits in the expansion over infinitely many orbits lfll l36ll37ll . 

Incorporating the results of the geometric ray-splitting orbits gives the following 
factor of the ray-dynamical approximation of the cluster determinant: 

The sum over r counts the repeats of the primary periodic orbits, the prime cycles p. If 
the logarithmic derivative with respect to ft) is taken, a result very similar to the spectral 
density for the interior problem is obtained. This is in agreement with the general result 
for the density of states for ray-splitting systems described in 113 Sh. Similar results for 
the case of flexural vibrations in the interior case are given in [39]. As the orbits are 
unstable and ip is symplectic, it is possible to expand, for each orbit, the instability 
denominator in (OTl in the inverse of its leading eigenvalue 1/Ap and to obtain a so- 
called Gutzwiller-Voros resummed zeta function similar to those of two-dimensional 
Hamiltonian flows lU]: 

oo 

Fg{z) = Y[;,\z). (22) 

where 

l/a(^) = n(l-4^) with 4^ = «, (23) 
p yi^pi^p 

8. SUMMARY 

Detailed studies of Helmholtz scattering determinants at small wave lengths have shown 
the influence of periodic orbits. The case of scattering from voids in two-dimensional 
elastodynamics was considered here with a discussion of the analytical contribution of 
periodic ray-splitting orbits to the scattering determinant. 
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